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We address the problem of pion production in low energy 7r-nucleus collisions. For the production 
mechanism we assume a simple model consisting of a coherent sum of single pion exchange and the 
excitation — followed by the decay into two pions and a nucleon — of the A r *(1440) resonance. The 
production amplitude is modified by the final state interaction between the pions calculated using 
the chirally improved Jiilich meson exchange model including the polarization of the nuclear medium 
by the pions. The model reproduces well the experimentally observed -n A — > uttA cross sections, 
especially the enhancement with increasing A of the tt + -k~ mass distribution in the threshold region. 



The past ten years have witnessed a considerable in- 
crease in our understanding of the tttt interaction. The 
success of chiral perturbation theory in describing 
the near-threshold behavior of the 7T7t amplitude is a 
clear success for the effective field theory approach to 
the problem. Unfortunately, at energies above a few hun- 
dred MeV, where effects of unitarity become important, 
chiral perturbation theory becomes unwieldy and other 
approaches must be used to obtain a good description of 
the free tttt scattering. One such approach is the meson 
exchange model developed by the Jiilich group |^] , which 
gives an excellent quantitative description of 7T7r and ttK 
scattering phases up to about f .5 GeV total cm energy. 
Another is the inverse amplitude method — a variant of 
the if-matrix approach — of Oset et al. |Q. Here we will 
consider the application of the former to the problem of 
pion production by pions on nuclei in order to investigate 
some of the predictions of the model for the behavior of 
the 7T7T scattering amplitude in the presence of a nuclear 
medium. The latest version of this model, which we refer 
to as the chirally improved Jiilich model || , respects the 
constraints on the 5-wave scattering length imposed by 
chiral symmetry while maintaining the quality of fit to 
the free tttt scattering data. We briefly summarize here 
the main features of the model and present some of the 
results of the model for tttt interactions in nuclear matter. 
For details of the model we direct the reader to ref. ||U ■ 

In the chirally improved Jiilich model, the interaction 
is driven by the exchange of p mesons plus contact in- 
teractions, as in the Weinberg lagrangian 0. This la- 
grangian is chirally symmetric in the massless pion limit. 
The Born approximation for this interaction is used as 
the potential in a three-dimensional scattering equation 
of the Blankenbecler-Sugar form |§. The solution of 
the scattering equation destroys chiral symmetry through 
both the partial summation of diagrams and the use of 
form factors, which are needed to ensure convergence of 



the integral equation. However, the off-shell behavior of 
the potential is prescribed in such a way as to preserve 
the scattering length constraint required by chiral sym- 
metry. 

The motivation to "chirally improve" the original 
Jiilich model resulted from studies of the behavior of 
the 7T7T interaction in nuclear matter as a function of 
density M. In these studies the polarization of the 
medium by the pion, principally through the production 
of nucleon-hole and A-hole configurations, significantly 
reshaped the 7T7t scattering amplitude. It was observed 
that models in which the potential term is not chirally 
symmetric in the limit of zero pion mass led to a rapid 
increase in attraction between the pions as a function of 
density, resulting ultimately in spontaneous 5*-wave pion 
pair condensation at densities about that of normal nu- 
clear matter. Further studies |I(||5|,^| showed that the 
onset of this instability could be shifted to higher densi- 
ties if the tttt potential term were required to satisfy the 
low-energy chiral constraints, or if the amplitude calcu- 
lated from the scattering equation could be forced to obey 
the same constraints. However, the build up of attrac- 
tive strength in the threshold and sub-threshold region 
persisted. In fig. |I| we show the 5-wave tttt scattering 
amplitudes at low energies in free space and at half nu- 
clear matter density using the chirally improved Jiilich 
model from ref. which will be employed in the fol- 
lowing. Whereas the isospin 1 = amplitudes (upper 
panel) undergo substantial reshaping in the medium, the 
1 = 2 channel (lower panel) is hardly affected. In the 
1 = channel, the combination of the attractive tttt in- 
teraction and the softening of the in-medium single-pion 
dispersion relation leads to a substantial accumulation of 
strength in the imaginary part of the in-medium ampli- 
tude (notice, however, that the real part is actually only 
enhanced for M TI < 300 MeV). In contrast, the purely 
repulsive interaction in the 1 = 2 channel prevents essen- 
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FIG. 1. Dimensionless S-Wave -kit scattering amplitudes 
for total isospin (upper panel) and 2 (lower panel) ; the dot- 
ted and dashed-dotted lines correspond to the real and imag- 
inary parts, respectively, in free space, and the long-dashed 
and full lines correspond to the real and imaginary parts, re- 
spectively, at a nuclear density of p — 0.5po- 

Since the in-medium scalar-isoscalar tttt amplitude is 
the main mediator of the intermediate range attraction in 
the nucleon-nucleon interaction, it is of prime importance 
for e.g. nuclear saturation. In fact, in refs. [ fH]]l2"| the 
use of the chirally improved Jiilich model in the Bonn 
Potential for the NN interaction [|l3| has been shown 
to be compatible with the empirical nuclear saturation 
point once medium modifications in the vector mesons 
exchanges (providing short-range repulsion) are also in- 
cluded. Here we are interested in a more direct assess- 
ment of the in-medium 7T7t S'-wave correlations. 

A challenging testbed for the in-medium modifications 
of the two-pion interaction discussed above has been set 
by recent experimental data taken by the CHAOS col- 
laboration @,||| at TRIUMF. An incoming pion beam 
of nominal kinetic energy T„.=282.7 MeV was directed 
at various nuclear targets, and two outgoing pions were 
detected in coincidence. With increasing nuclear mass 
number, the corresponding invariant mass spectra show 
a strong enhancement of 7r + 7r~ pairs just above the two- 
pion threshold of M W7r = 2m T , whereas only very minor 
variations are observed for 7r + 7r + pairs. 

That the 7r + 7r~ strength in the scalar-isoscalar chan- 
nel in a nuclear medium can be strongly enhanced at 
threshold was predicted in ref. fig] , based on inclusive 
data from the CHAOS collaboration [fL7|. A first indica- 



tion that this threshold enhancement can account for the 
measured (7r + , 7t + 7t~) cross section at low invariant mass 
was obtained from a schematic calculation in ref. [p| . A 
clue to the possible underlying mechanism for this dra- 
matic effect is provided by the isospin decomposition of 
the charged states. For S-waves one has 



Mn+ir- = -M 



1=0 
1=2 



(1) 



showing that while the 7r + 7r + channel is pure isotensor, 
the 7r + 7r _ channel is predominantly isoscalar. Thus, the 
in-medium 7T7t amplitudes § shown in fig. [l] clearly ex- 
hibit desirable features: little modification of the scalar- 
isotensor amplitude, and a significant low energy en- 
hancement of the scalar-isoscalar amplitude. To quan- 
tify the comparison, the actual cross sections for the 
ttA — > irnA production process must be calculated, in- 
cluding the experimental acceptance of the CHAOS spec- 
trometer, which is crucial to reproducing the precise 
shapes of the observed spectra. 

Since our emphasis here is on possible medium effects 
in the -k-k interaction, we will treat the pion production 
process in a simplified manner. First, we assume that it 
always proceeds as an elementary process on a single nu- 
cleon. Second, we account only for the two most impor- 
tant contributions which, according to refs. Jlj|[l5|], are 
the one-pion exchange (OPE) reaction (contributing to 
both isoscalar and isotensor channel) and the TV* (1440) 
resonance formation (leading to isoscalar 7T7t states only) . 

Let us first discuss the OPE contribution. The corre- 
sponding 4-diffcrential cross section is given by 
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where M nv is the two-pion invariant mass, D v (t) 



_1 the propagator of the exchanged pion, pi a f, the 



incoming pion laboratory momentum, and q, q' , 
the in/outgoing pion momentum and their relative angles 
in the two-pion cms (/ w n n =1.01 and IF =2 for charged 
pion exchange). Azimuthal symmetry in the 7T7t inter- 
action makes the <p cm integration trivial. Furthermore, 
since our focus will be on the near threshold region, 
the iS-Wave contribution is expected to be dominant. 
This is supported by experimental angular distributions 
in cos# cm which do not show any trace of other than 
(isotropic) S'-wave components [|l4j . We will thus neglect 
P- and higher partial waves. The invariant mass spectra 
are then given by 
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For a meaningful comparison with the experimental 
spectra it is essential to include the acceptance of the 
CHAOS spectrometer, represented by an 'acceptance fac- 
tor' Acc(t, M^ir, Stot) in eq. (||). We evaluate it by Monte 
Carlo techniques: in the two-pion CMS we randomly gen- 
erate pairs of angles (</> cm , cos(9 cm ), uniformly distributed 
over the sphere due to the 5-wave nature of the tttt- 
amplitude. (To include higher partial waves one would 
need to weight the polar angles with the corresponding 
angular distribution obtained from the total amplitude 
for the particular M wv ,q,q'.) For given t, M^, s to t one 
can then determine the total momentum P of the pion 
pair in the laboratory frame and apply the corresponding 
Lorentz boost to transform the two generated tracks into 
the lab system. In the latter, the experimental accep- 
tance cuts can then be readily applied, i.e. — 0° ± 7 
or = 180° ± 7 for the out-of-reaction plane opening 
angle of the pion pair and T ± = 10° — 170° for the in- 
plane angle of both single pion tracks accounting for the 
dead regions around the beam-direction. At fixed t, M W7r 
and Stot the acceptance probability is then determined 
as Acc = N acc /N to t, with N to t the number of trials and 
N aC c the number of events that fall into the experimental 
acceptance. 

In order to assess the A* (1440) contribution we model 
the two-pion production amplitude as proceeding via 
an intermediate scalar-isoscalar resonance a of invari- 
ant mass -Mn-Tr, which subsequently decays into two pi- 
ons, i.e ttN -> TV* (1440) -> dN -> ttttN. Since in the 
Jiilich model there is no low-lying genuine a resonance 
(the 'cr(550)' being chiefly generated through attractive 
i-channel p exchange, and the /o(980) being realized as a 
KK bound state) , we identify the a with the e resonance 
located at around M W7r ~1400 MeV. Its coupling to tttt 
states is well determined by a satisfactory fit to the 
phase shifts above 1 GeV |§. For the A*(1440) Ae vertex 
we assume scalar coupling, 



N* Ne = gN*Ne *jv* 6 * 



- N 



h.c. 



(4) 



and estimate the corresponding coupling constant g^'Ne 
from the experimentally measured branching ratio of 5- 
10% for A* (1440) -> N(wir)jj% ave @. Similarly, for 
the entrance channel, we employ the usual interaction 
vertex 



J~N*Ntt — N* 



3 ■ q T a TT a * w + h.c. , (5) 



and adjust the N*Ntt coupling constant to the branching 
ratio for A* (1440) -> Avr of 60-70% @. In Born ap- 
proximation, the cross section for ttN — ► A* (1440) — > eA 
with an e of invariant mass M T „ and 3-momcntum P is 
then obtained as 



&TrN^N'^eN(stot, M-xt:) 
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with the spin averaged invariant matrix element in the 
ttN cms 

IXi 12 — ^ fN'Nii 2 ft" 
\MN*eF ~ -7T n g N *Ne 



(7) 

The decay of the e is now included by folding the cross 
section eq. ([j]) over the fully dressed e spectral function 
which automatically accounts for the final state interac- 
tion of the pion pair jiy] : 

VirN^N*^tnrN(stot) = 

— — A e (M 7rff ) (TnN->N*->eN(Btati M^) ■ (8) 



with 

MM**) = (D° £ (M^)f tt 4 

x [1 -7T q' M wn \mMZ{M^)] , (9) 

where v tlT7T denotes the ctttt vertex function and the 
full tttt scattering amplitude in the JI = 00 channel. Af- 
ter a variable transformation, dP = E e (P)/(2MNP)dt, 
and inclusion of the experimental acceptance, we finally 
obtain the differential mass spectrum 
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The cross-section we wish to calculate requires a co- 
herent sum of the OPE and A* (1440) amplitudes. The 
result of this yields the M T?r differential cross section 
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where the amplitudes M° PE and Ai N can be read off 
by comparison with eqs. (||) and (|l0|), respectively. 

In the upper panels of figs. ^ and || our results for 
deuterium target are displayed. The cross sections are 
calculated in absolute units using the free tttt interac- 
tion. The full lines correspond to a calculation with the 
Monte Carlo-simulated acceptance as described above. 
In the tt + tt~ channel the theoretical curves somewhat 
overestimate the data for M vv < 320 MeV, and under- 
estimate them for > 350 MeV, a feature that is 
also shared by the model calculations shown in ref . jl5| ; 
in the tt + tt + channel the agreement is quite satisfactory. 
To check our determination of the experimental accep- 
tance we also performed calculations with an 'effective 
acceptance' (dashed lines in the upper panels); this was 
extracted from a fit to a phase space simulation given in 
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ref. which, after dividing out the phase space, leads 
to a simple M^-dependent acceptance factor Acc{M vv ). 
Whereas the description of the tt + it~ channel turns out 
to be somewhat improved, the opposite trend is found 
in the 7r + 7r + channel. In general, both ways of imple- 
menting the experimental acceptance agree reasonably 
well — within 20 percent. 

When moving to finite nuclei, the effect of the nuclear 
Fermi motion has to be included. In principle this is very 
complicated, even if everything about the location of the 
interaction within the nucleus were known. We adopt 
instead a simplified procedure that, as we will see, can 
adequately account for the observed cross-sections. First, 
we assume that the production of the pion pairs occurs 
at a definite density. In principle one should then aver- 
age over all momenta within the Fermi sphere for that 
density, which would mean calculating the cross-section 
for non-collinear collisions, and we should also account 
for the effective mass of the nucleon as a function of mo- 
mentum. This is computationally tedious and, given the 
precision of the data, largely unnecessary. One can show 
that sufficient accuracy — on the order of one percent or 
better — is obtained if the average is performed only over 
the component of the nucleon momentum along the beam 
direction. In this case one has 
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with s tot (k z ) = ml + Mfj + 2(u piab E N (k z ) - pi ab k z ) and 
the Fermi momentum given in terms of the local den- 
sity by pn = 2kp/3n 2 . (In fact, we adopted the same 
procedure for the deuteron calculations with an 'effec- 
tive' Fermi momentum representing the kinematical limit 
where the two nucleons recoil together.) 

Of course pions undergo substantial absorption in the 
nuclear medium — an effect which, as implied in the pre- 
vious paragraph, we did not include in the present analy- 
sis. Here we address instead the question of whether the 
medium modifications of the final state 7T7t interaction 
can account for the dramatic reshaping of the mass dis- 
tribution observed in the CHAOS experiment. To fix the 
units for the in-medium results we assume that our the- 
oretical total cross sections account for the same fraction 
of the measured cross sections as on the 2 H target (note 
again that our deuteron results are in absolute units), 
which ranges between 80% and ~100%. The use of the 
in-medium rnr amplitudes in the differential cross sec- 
tions then leads to the results shown in the lower panel 
of figs. ^| and |^, where we have employed the parameter- 
ized, 'effective' experimental acceptance. 

In the 7r + 7r + channel the density dependence is very 
weak, the major effect being a smearing of the cross sec- 
tion due to the nuclear Fermi motion. However the tt + tt~ 
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FIG. 2. Results for pion production cross sections in the 
■r + tx" channel using the chirally improved Jiilich tvtt interac- 
tion from ref. J5I; upper panel: on the deuteron target using 
the Monte Carlo-simulated acceptance (full lines) and the pa- 
rameterized, 'effective' acceptance; lower panel: on Calcium, 
using the effective acceptance, for two different densities. The 
2 H and 40 Ca data are from ref. Jl^| and respectively. 
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FIG. 3. Results for pion production cross sections in the 
7r + 7r + channel using the chirally improved Jiilich -kit interac- 
tion; identical line identification as in fig. ^| 

channel exhibits very substantial medium effects; in par- 
ticular, the strong threshold enhancement observed in 
the data is nicely accounted for by our model, being en- 
tirely due to the effects in the in-medium isoscalar iS-wave 
■kit amplitude . The average nuclear density which 
seems to reproduce the 40 Ca data best turns out to be 
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slightly above half the saturation density (p = 0.5-0.7po)- 
This is in line with the average density of 40 Ca and the 
fact that the two pions are probably created near the 
surface. One should, however, realize that the two out- 
going pions are at rather low energies and can therefore 
pass the nucleus with little attenuation. Also note that 
the ratio of the threshold peak to the main maximum at 
A/tttt ~ 350 MeV increases within rather small increases 
of density, a trend that is also clearly seen in the exper- 
imental data when going to heavier nuclei (from 12 C to 
40 C a to 208 Pb) p5| , where one would expect to probe 
slightly increasing average densities. 

To summarize, we have analyzed A(tt + , n + n^) pro- 
duction cross sections on the deuteron and nuclei at 
low energies using a realistic model of the tttt interac- 
tion which accurately describes free rnr scattering over 
a wide range of energies and incorporates 'minimal chi- 
ral constraints'. Our emphasis was on the S-wstve irir 
final state interaction and, in particular, its modifi- 
cations by the medium. The latter were included in 
terms of the standard renormalization of pion propa- 
gation through nucleon/A-hole excitations entering a 
Lippmann-Schwinger-type equation for the tttt ampli- 
tude. Many aspects of the complicated reaction dynamics 
have been treated approximately — e.g. the pion produc- 
tion process, the nuclear Fermi motion and absorption — 
or neglected altogether — e.g. Coulomb interactions and 
the dependence of the medium corrections on the total 
two-pion three-momentum. Nevertheless we find that our 
calculations, including experimental acceptance, agree 
reasonably well with recent data taken by the CHAOS 
collaboration at TRIUMF in both the ir + ir~ and the 
7r + 7r + channel. In particular, the strong enhancement 
of the 7r + 7r~ cross section very close to the free two-pion 
threshold can be well reproduced by ascribing it to the 
medium-modified scalar-isoscalar tttt interaction. Such 
an enhancement is absent in the isotensor ir + ir + interac- 
tion, which is also consistent with the data. The physical 
origin of these features is rather simple: the softening of 
the in-medium single-pion dispersion relation acts coher- 
ently with the attractive interaction in the isoscalar chan- 
nel, but is not effective in the purely repulsive isotensor 
channel. 

The outgoing two pions being on-mass-shell precludes 
any statement as to whether the threshold enhancement 
of the in-medium isoscalar tttt amplitude entails substan- 
tial subthreshold effects at higher densities, as predicted 
in various tttt models. However, since the Tr + n~ correla- 
tions in the J—I—0 channel are of great importance for 
the understanding of nuclear binding, further experimen- 
tal and theoretical studies of the subject will certainly be 
of great interest. 
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